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Abstract 

In this paper we consider a resonance problem, in a generic regime, in the con¬ 
sideration of relaxation of ground states of semilinear Schrodinger equations. 
Different from previous results, our consideration includes the presence of reso¬ 
nance, resulted by overlaps of frequencies of different states. All the known key 
results, proved under non-resonance conditions, have been recovered uniformly. 
These are achieved by better understandings of normal form transformation and 
Fermi Golden rule. Especially, we find that if certain denominators are zeros (or 
small), resulted by the presence of resonances (or close to it), then cancellations 
between terms make the corresponding numerators proportionally small. 
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1 Introduction 

We consider the following 3-dimensional semilinear Schrodinger equations 

idti/)(x,t) = - A 4>{x,t) +V(x)il>(x,t) + \ip{x,t)\ 2 ip(x,t), (1) 

il>(x, 0 ) =ipo{x) € 7t 2 (M 3 ) small. 

where V : ffi 3 —> R is the external potential. 

Such equations arise in the theory of Bose-Einstein condensation, nonlinear 
optics, theory of water waves and in other areas. 

We start with formulating the problem. The potential V : K 3 —> R is a 
smooth, and rapidly decaying function. And if it is trapping potential, namely 

inf (/, (-A + F)/) = -e 0 <0, ( 2 ) 

Il/ll2=i 
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then this linear unbounded self-adjoint operator — A + V, mapping C 2 into £ , 
has a ground state (j> € C 2 with eigenvalue —eo- Moreover the eigenvalue must 
be simple. Besides the ground states the linear operator might have some other 
finitely many neutral modes with nonpositive eigenvalues — e*,, k = 1, • ■ ■ , N. 
Its continuous spectrum spans the interval [0, oo). It is well known that for the 
type of potential V we chose, there is no positive eigenvalues, see e.g. Ba¬ 
in the nonlinear setting, the ground state bifurcates into a family of solitary 
wave solutions, see e.g. [21], 

ip(x, t) = e lXt (j) X (x) (3) 

with AeM being close to eo and </> A = C^J |eo — A| <j> + O(\eo — A| ®). 

There is a rich literature on studying the orbital stability and asymptotic 
stability of the soliton manifold. By results in [251 [26] B3l it is well known that 
that the ground state manifold is orbital stable in the H 1 space. After these, 
many attempts were made on proving the asymptotic stability of the ground 

state manifold, see e.g. 0 iisiisi 123100 in EDi Ezmanai. 

In [5] El 0 [24j [23j [20] , it is assumed that the linear operator — A + V has 
a ground state (with eigenvalue —eo < 0), and only one simple neutral mode 
with eigenvalue — e\ < 0 satisfying 2ei < eo- In [221 . multiple neutral modes was 
considered. Their eigenvalues — e*,, k = 1, ■ • • , N must satisfy two conditions: 
(1) 2efc < eo, and (2) the so-called non-resonance condition, namely there do 
not exist rik £ Z, k = 0,1, • • • , N such that 

H l nfc l ^ 0 and ^ rikek = 0, (4) 

k 

see also the non-resonance conditions for multiple neutral modes in Id. 

On the technical level, the condition ([4]) was needed to prevent small de¬ 
nominator from appearing. 

In [HI E], the author, together with M. Weinstein, improved the above 
results by studying degenerate neutral modes, i.e. 

efe = e x , fc = 1,2, ■ ■ ■ ,N, (5) 


or nearly degenerate. 

The main purpose of the present paper is to include the presence of the 
resonance, specifically by removing the conditions ([4]) and (0, and to show all 
the proved results still hold, uniformly. 

A graphic illustration is in Figure [I] 

On the technical level, we achieve this by re-defining normal form transfor¬ 
mation and Fermi Golden rule. Especially we show that if some denominators 
are small, caused by the resonances, then their corresponding numerators are 
proportionally small, resulted by cancellations between terms. For the details, 
we refer to Sections 0 and El and Proposition 111.11 

To the best of our knowledge, our result and techniques are new. 

Related problems. The main motivation of the present work is to under¬ 
stand certain types of small divisor problems. In many literatures, for example 
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Figure 1: spectrum of — A + V 


in in studying problems arising from statistical mechanics and PDE, 

the small divisor was avoided, by choosing initial conditions, for example. 

The technical advantage in studying the present problem is that we only need 
to expand the solutions, in certain small parameters, finitely many times (two 
times), instead of infinitely many times as in [3]. We expect the normal form 
transformation invented here, which makes it easy to see the crucial cancellations 
between terms, together with algebraic structures observed in [9U271 m, for 
higher order iterations, can be applied to the other problems. These will be 
addressed in subsequent papers. 

Structure of the paper. The paper is organized as follows. Some basic 
properties of the equation are studied in Section [3] The linearized operator 
L( A), obtained by linearizing the ground state solution, is studied in Section [2 
The Fermi Golden rule condition is in Section [5] The main Theorem is stated 
in Section [6] The decomposition of the solutions, and the governing equations 
for various parameters and functions are in Section 0 The main Theorem is 
reformulated into Theorem 18.11 in Section [8] The proof of different parts of 
Theorem 18.11 are in Sections M Col and cu The main theorem is proved in 
Section C2 

Most parts of the paper follow the steps in the previous papers mm- 
Technically the main differences are in normal form transformation in Section 
[HU and subsequent sections. 
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2 Notation 

(1) a+= max{a, 0}, [r] = max? ez {f < t} 

(2) 5ftz = real part of z, 3z = imaginary part of z 

(3) Multi-indices 


z = (zi,...,z N ) S C N , z = (zi,... ,z N ) (6) 

a € = zl 1 ■■■z a N N 

l a l = l a il + • • • + M 

(4) H s = Sobolev space of order s 

(5) 



(6) aess{L) = a c (L) is the essential (continuous) spectrum of L , 
ad = C — <J C {L) is the discrete spectrum of L. 

(7) Pd(L) bi-orthogonal projection onto the discrete spectral part of L 

(8) P C {L) = I — Pd{L ), bi-projection onto the continuous spectral part of L 

(9) (f,g) = f f 9 

(10) ||/||p = L p norm, 1 < p < oo. 
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3 Basic Properties 

Equation m is a Hamiltonian system on Sobolev space 'H 1 (K 3 ,C) viewed as a 
real space ’H 1 (R 3 ,R) © 7d(R. 3 ,K). The Hamiltonian functional is: 

m ~ J\\{ w\ 2 +v\ipf)+^\% 

Equation m has the time-translational and gauge symmetries which imply 
the following conservation laws: for any t > 0, we have 

(CE) Conservation of energy: = £(ip(0)); 

(CP) Conservation of particle number: 

= Af(4>( 0)), where A f(ip) := f \fi\ 2 . 

In what follows we review the results of the existence of soliton and their 
properties. 

The following arguments are almost identical to those in mmm except 
that here we have multiple neutral modes (or excited states), hence we state 
the results without proof. We assume that the linear operator — A + V has the 
following properties 

(NL) The linear operator —A + V has eigenvalues — eo < —ei < • • • < —ejv 
satisfying eo < 2ei. —eo is the lowest eigenvalue with ground state <f> > 
0 , the eigenvalue — ei, ••• , — ejv might be degenerate with eigenvectors 

clin clin clin 

SI 5 S2 j * * *> S N m 

In the nonlinear setting the ground state bifurcates into a family of solitary 
wave solutions of o , see e.g. [24] . 

= e lXt cl) X (x) (7) 

and the function (f> x > 0 has the following properties, see e.g. [241 . 

Lemma 3.1. Suppose that the linear operator — A + H satisfies the conditions in 
[NL] above. Then there exists a constant 6 q > 0 such that for any A £ [eo — <5, eo) 
0 has solutions of the form ) = e lXt (f x {x ) £ C 2 with 

<t> x = 5<j> + 0(<5 3 ) (8) 

and 5 = (f cf> 4 (x)dx) 2 (e 0 - A) 1 / 2 + o((eo - A) 1 / 2 ), moreover 
d\<j) x = O((e 0 - A) -1/2 )</> + o((e 0 - A) 1/2 ). 

4 The Linearized Operator 

After linearizing the solution around solitary wave solution 0, namely consid¬ 
ering the solution of 0 ip(x,t) = e lXt [(j) X {x) + R(x,t)\, then the linear part of 
the equation for R(x, t) is 

d t R = L(X)R (9) 
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with R \= ( ReR, ImR) T , and the linearized operator L{ A) is defined as 


L( A):= 


0 L_( A) 

-i+(A) 0 


(10) 


with L±(A) being linear Schrodinger operators defined as 

L-(X) := -A + V + A + (<^ A ) 2 , and L+(A) := -A + V + A + 3(<^ A ) 2 . 

By general result (Weyl’s Theorem) on stability of the essential spectrum 
for localized perturbations of J(—A) [T7j . 

cr ess (T(A)) = (— ioo,— zA] D [*A, ioo) 

if the potential V in Equation ([1]) decays at oo sufficiently rapidly. 

Next we study the eigenvalues and eigenvectors of L{ A). The proof can be 
found in m, hence we omit it. 

Lemma 4.1. Let L{ A), or more explicitly, L(X(S),6) denote the linearized op¬ 
erator about the the bifurcating state (/> A , A = A(<5). Note that A(0) = eg. 


It has an eigenvector 
eigenvalue 0: 


L{ X) 


0 


0 


= 0 , 


and an associated eigenvector 


d\4>* 

0 


L( A) 


d\ <t> x 
0 


Corresponding to the (possibly degenerate) eigenvalue, — e\, —e 2 , 
of — A + V, the matrix operator 


with 


( 11 ) 


-e-N, 


L(A = eg, 6 = 0) 


has eigenvalues 


±iE n (e 0 ) = ±i(e 0 - e n ), n = 1,2, • ■ • ,N. (12) 


For S > 0 and small, these bifurcate to (possibly degenerate) eigenvalues, of the 
operator L(X), ±iE\(X),..., ±iEn(X) with eigenvectors 


( ±i m ) ' ( ±% 2 ±iVN ) 

with t( n , rj n being real functions and 

‘ £n) — i ^?n) — 0, (fni Pro) — 


(13) 


(14) 


Moreover, for S sufficiently small 2E n (X) > A, n = 1, 2, • • •, N, (resonance at 
second order with radiation). 

Furthermore we need the following condition on the threshold resonances. 
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Definition 1. A function h is called a threshold resonance function of — A + V 
at 0, the endpoint of the essential spectrum, |h(x)| < c(x) _1 + and h is C 2 and 
solves the equation 

(—A + V)h = 0. 

A function h is called a threshold resonance function of L( A) at p = ±*A, 
the endpoint of the essential spectrum, \h(x)\ < c(x}“ 1+ and h is C 2 and solves 
the equation 

(L( A) — p)h = 0. 

In this paper we make the following assumption: 

(SA) —A + V has no threshold resonance at 0. 

This assumption is generic since it is known that the threshold resonance is 
unstable, see e.g. m- Based on this assumption it is well known that 

Lemma 4.2. If |eo — A| is sufficiently small and the assumption (SA) holds, 
then L(X) has no threshold resonances at p = ±*A, and L( A) has no other 
eigenvectors and eigenvalues besides the ones listed in Lemma \4.1\ 

We denote the projection onto the essential spectrum of linear operator L( A) 
is Pf = 1 — P£. In the following we give the explicit form of the projection Pd, 
whose proof for N = 1 can be found in m, the proof of the general cases are 
similar, hence omitted. 


Proposition 4.3. For the non self-adjoint operator L(X) the (Riesz) projection 
onto the discrete spectrum subspace of L(X), Pd = Pd(L(X)) = P(f, is given by 


Pd 


cAIIPII 2 



N 




tin 

illn 


-ir]n 

tin 


f,n \ / iVn 

iVn / \ fn 


(15) 


We define the projection onto the continuous spectral subspace of L( A) by 


P c = P C (L(X)) = P C A = I - P d . (16) 


5 The Negativity of Fermi Golden Rule in Ma¬ 
trix Form 

The Fermi Golden rule plays an essential role in determining the decay rate of 
the neutral modes. For simple neutral modes, as in [2U EH 0 m, the form is 
simple because the number is only one, or if more than one then they can be 
separated into independent ones by a near identity transformation (see )22l f8]). 
The problem of multiple neutral modes is more involved due to the fact that 
multiple coupled parameters appear and they can not be separated. 
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Next we define the new Fermi Golden Rule condition. 

Define a function e : t 3 x R 3 -> € as 

e(x,k) := [l + (-A x -\k\ 2 -iO)~ 1 V(x)\ 1 e lk ' x . (17) 

It is known that for any fixed k ^ 0, this function is well defined by the type of 
potential V we chose, see e.g. m- And it satisfies the equation 

[~A X + V(x) - \k\ 2 ] e(x, k ) = 0. 

We define complex functions 5' m n on the 2-dimensional unit sphere § 2 . For 
any a G S 2 , 

^m,n(cr) = 4' n , m (cr) := f e(x, \k\ m . n a)(/)(x) (x)£, 1 ™ (x) dx (18) 

It 3 

with G R + defined as 

:= (co C m e n ) 2 , 

where, recall that we assume that 2e t < eo, l = 1,2,- ,N, and </>, fffff are 
eigenvectors of —A + V with eigenvalues —eo and —e m , m = 1, • • • , N. 

Now we state our Fermi-Golden-rule condition. 

(FGR) For any scalar vector 3 = (z\, Z 2 , • • • ,z n ) € C N satisfying \z\ = 1, the func- 

N 

tions defined on the unit sphere § 2 , ’I'm,ri(c r )z m z„ is not identically 

m,n—l 

zero. 

Its important ramification is that there exists some constant C > 0, such 
that for any z = (z i, z 2 , - - - , 2 n ) G C N , 


r(z, z) := || ^2 ^m,n(^)ZmZn\\c 2 (S^ > C|z| 4 . (19) 

m, n 

Here we use that Ev n{&)zmZn is smooth in <r, hence if it is not identically 

m, n 

zero, we have the estimate above. 

Remark 1. If the set of eigenvalues {e k | k = l,--- , N} can be grouped into 
well separated clusters, namely 

{e k \ k = !,-■■ ,N}=U l A l (20) 

with properties that, for some constant Cq = 0(1), 

|e„ - e m | > c 0 if e n G A u e m G A k with l £ k. 

Then Fermi Golden rule assumption can be relaxed: for each fixed l, the each of 
the functions If ’m,n(o') z m z n is not identically zero, with £ M 2 X 
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6 Main Theorem 


In this section we state precisely the main theorem of this paper. The key fact 
is that despite of the possible presence of resonance, the main results remain 
the same as in mm- 

For technical reasons we impose the following conditions on the external 
potential V of ©: 

(VA) V decays exponentially fast at oo. 

Recall the notations £ = (£ 1 , • • • , £jv) and r) = (rj 1 , • • ■, ? 7 at) for components 
of the neutrally stable modes of frequencies ±iE n (\), n = 1 ,■■■ ,N, of the 
linearized operator L( A) defined in (fTOl) . 

Theorem 6.1. Suppose that Conditions (NL) in Section [H (SA) in Section 
[U (FGR) in Section 0 and (VA) above are satisfied. Let v > 0 be fixed and 
sufficiently large. 

Then there exist constants c,e o > 0 such that, if 


inf 

7£R 


V-o - e* 7 



(Re z^) ■ f + i (Im z^) ■ r^j 


'H 3 ’‘ 


then there exist smooth functions 


< 


c | 2 (0) | < eo, 


( 21 ) 


A (t) : K + -> X, 7 (i) : M + ->• R, z(t) : R+ C d , 
R(x,t) : M. d x R + ->• C 


such that the solution of NLS evolves in the form: 
fi>(x,t) = e if o Hs)ds e i 7 ( t ) 

x [cj) X + ai(z, z)d\4> x + ia 2 (z, z)fi x + (Re z) ■ £ + i(Imz) ■ i] + R], 

( 22 ) 


where limt_ i . 00 A (t) = Aoo, for some Aoo £ T. 

Here, a±(z, z), a 2 (z,z ) : C N x C N —¥ R. and z — z : <C N x — > C N are polyno¬ 
mials of z and z, beginning with terms of order \z\ 2 . Their explicit definitions 
will be given in (1581) . 

Moreover: 

(A) |z(t)| < c(l +t)~i and, there exists a polynomial F(z,z) = 0(|z| 4 ) £ R 
such that z satisfies the initial value problem 

d t [\z\ 2 + F(z, z)} = -CT(z, z) + 0(( 1 + 1)-%) (23) 

where C > 0 is a constant, T(z,z) = 0(\z\ A ) is a positive quantity defined 
in (fl9l) . 
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(B) R(t ) = (ReR(t), ImR(t)) T lies in the essential spectral part of L(X(t)). 
Equivalently, R(-,t) satisfies the symplectic orthogonality conditions: 

u:(R,i(j) X ) = uj(R,d\(j) X ) = 0 

iv(R, iijn) = u(R, £„) =0, n = 1,2, • ■ N, (24) 

where uj(X,Y) = Im f XY and satisfies the decay estimate: 

IKl + sV^WIIa^ca + lil)- 1 . (25) 

The main theorem will be reformulated into Theorem 18. II below. 


7 The Effective Equations for z, A, 7 and R 


In this section we derive equations for i, A, 7 and R. 

We decompose the solution as 

ip(x,t) = x ( s ) ds e il,(t) 

N 

X (j) X + ai(px + + Pn)(,n + i (fin + 9n)^n + R 


N 


n= 1 n—1 

= e* to A ( s ) ds e *T(t) _|_ ai ^)A _|_ j a2 ^A _|_ ( a -|- p) . £ j(^ + g) ■ I? + i?] 

(26) 


Here and going forward, we’ll use the notations: 

(X (op, • ■ • , CK^v) ; (d (/A > ■ ■ • ) /Av) ; 

£ = ( 6 , • ■ • , 6 v) T , V = (Vi ,---,Vn) T - 


Introducing 


2 = a + ifi, 


we have 

a = 2 < ' 2 + ^’ ^ 

and we seek aj = aj(z,~z) = 0(\z\) 2 , j = 1,2 and pj = pj(z,z ) = 0(\z\ 2 ), 
polynomials in z and z, which are of degree larger than or equal to two, and 
are real. Substitution of the Ansatz (1261) into NLS, equation (fl|). we have the 
following system of equations for R , defined as 


R := 


) , R\ ■= ReR , R 2 := ImR : 
El ) 
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as 


t>,R=LmR + iJR-jm*)-( X+tal-h] ) 


( £, ■ [E{\){/3 + q) - d t {a +p)\ \ 

V -V ■ [E{X){a + p) + d t (/3 + q)\ ) 

, ( (P + q) ■ V \ _{( a\dl<i) x + (a + p) ■ d x £; \ 
V -(«+p)-£ ) V + (P + q) ■ d x r] ) 


(27) 


Here, 


with 


JN(R, z) := 


ImN(R,z) \ 
-ReN(R,z) ) 


(28) 


ImN(R, z) := |<^ + h + ihfh - ((/> X fh 
ReN(R, z ) := [|<^ A + h + ih? - (</> A ) 2 ] W> A + h) ~ 2(</» A ) 2 J 1 
I\ = A± + A2 + i?l, I2 = B\ + B 2 + i ?2 
-4i = a • £, ^2 = aid\4> x + p-£, 

Bi = P ■ rj, B 2 = a 2 (t> X + Q'V- (29) 

From the system of equations (l27l) and the orthogonality conditions (fT4l) and 
(l24l) we obtain equations for A, 7 and z n = a n + ip n , n = 1,..., N : 


d t (a n +p n ) - E n (X)(P n + q n ) + (ImN(R, z), rj n ) =F ln ; 

dt(Pn + q n ) +E n (X)(a n +p n ) - (ReN(R,z),£ n ) =F 2n ; 

1 


7 + d t a 2 - on - ' . ( ReN(R , z), (j) x ) =F 3 ; 

XY 1<P a ) 

1 


A + dt<i\ + 


(<A A )^a) 


(ImN(R, z), 4> x ) =F a . 


(30) 

(31) 

(32) 

(33) 
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Finally, the scalar functions F hn . j = l,2,Fs,F4, are defined as 

Fin :=7 {(P + O) ■ V,Vn) - \ai{d 2 x (j) x ,r] n ) - A {(a+p) ■ d x ^rj n ) 

- j{R 2 ,ri n ) + \(Ri,d\r) n ), 

F 2n ■= - 7 {(a+p) ■ £,£„) - Aa 2 {(/>x,£n) - A((/3 + q) ■ d x r/,£ n ) 

+ 7(-Rl, £n) + A(i?2, d X t; n ), 

1 

F, := 


F 4 := 


(0\0A> 

A(F 2 ,</>aa) - i( R u<t>\) - (7 (a+p) • £ + Aa 2 <^ + A(/3 + g) • d x r],<t> x ) 
1 


(0 a >0a> 


A(i?i,^) + 7 (F 2 , </> A ) + < 7(/3 + 9 ) • 7? - Xaid 2 cj> x - A(a + p) • 5 a£, <^ A ) 

(34) 

To facilitate later discussions we cast and (ED into a convenient form. 
Since a n and /3 n are real parameters, it is equivalent to study the complex 
parameters z n := a n + i/3 n . Compute (l301) +i (15Tl) to find 

d t z n + iE n (X)z n = - d t {pn + iq n ) - iE n (X)(p n + iq n ) - ( ImN(R 1 z), p n ) 

+ i(ReN(R, z), £„} + -fi„ + iF 2ri . (35) 

Note that (1501) and (1511) can be recovered from the equation above by considering 
the real and imaginary parts of (1351) . 

Remark 2. (a) Recall the estimate of Remainder in (l48l) . By (1501) - (l33l) we 

have 

A, 7 , dtz n + iE n (X)z n = 0(\z\ 2 ) + Remainder. (36) 

(b) The functions aj(z,~z),j = 1,2, p n (z,~z ), q n (z,~z),n = 1 will be 

chosen to eliminate “non-resonant” terms: z a ~z b , 2 < a + b < 3. 

Finally, we derive an equation for 

R = P x{t) R = P C R, 

the continuous spectral part of the solution, relative to the operator, F(A(t)). 
Applying P c = P X ^ to (l27l) to use that (see (flSll l 


±ir]n 


= Pr 


0 


= Pr 


d X <t> } 

0 


= 0 


to remove many terms on the right hand side, and using the commutator iden¬ 
tity: 


P c d t R = d t R - Xd x P c R 


(37) 
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we obtain 


d t R = L(X(t))R - PWJN(R, Z ) + L^R + g. (38) 

Here the operator L ^ ^ and the vector function Q are defined as 

L (jvV) = A9 A P c A W+7P c A(t) J, (39) 

G = pHt) ( + Xa } dl(j) x - A {a + p) ■ d x f \ 

° V -Hot + p) ■ £ - Xa 2 <fc - \(P + q) ■ d x r) )' 

8 Reformulation of The Main Theorem 

The proof of Theorem 16.11 we use the following result, which gives a more 
detailed characterization of the terms in the decomposition. 

Theorem 8.1. The function R in (1221) of Theorem 16. /I can be decomposed as 
R= J2 tf™,»(A )z m z n + R (41) 

|m| + |n|=2 


where R m ,n are functions of the form 

Rm,n = [L{\) + lE( A) • {m - n) - 0] Vm,n 

<t>m,n are smooth, spatially exponentially decaying functions. The function R 
satisfies the equation 


d t R = L(X)R + M 2 (z,z)R + P c N 2 (R,z) + P c S 2 (z,z), (42) 


where 

(1) S 2 (z,z) = 0(\z\ 3 ) is a polynomial in z and z with A -dependent coefficients, 
and each coefficient can be written as the sum of functions of the form 

[L( A) ± i(E m ( A) + E n { A)) - O]- fc P c 0 ±fc (A), (43) 

where k = 0,1, 2 and the functions <f>±k{ A) are smooth and decay exponen¬ 
tially fast at spatial oo; 

(2) M 2 (z,z) is an operator defined by 

M 2 (z , z) := fP c J + XPcX + X, (44) 

where X is a 2x2 matrix, satisfying the bound 

\X\ < c\z\e- e °W. 
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(3) N 2 (R, z) can be separated into localized term and nonlocal term 

N 2 = Loc + N onLoc 


(45) 


where Loc consists of terms spatially localized (exponentially) function of 
i£l 3 as a factor and satisfies the estimate 


IKzH-A + l)Loc\\ 2 + ||£oc||i + \\Loc\\ 4 < c(|^| 3 (i) + |z|(i)||(:r) V R\\ 2 ) 

and NonLoc is given by 

NonLoc := ( R\ + R\)JR. 

Here v is the same as in Theorem 16.11 


(46) 

(47) 


In the rest of the paper we denote by Remainder(t) any quantity which 
satisfies the estimate: 


\Remainder(t)\ < |z(i )| 4 + \\{x)~ v R{t)\\l + II^WIlL + \*V)\ U^mh- 
The functions A, 7 , z have the following properties 

A = Remainder (f); 


(48) 


(A) 

(B) 


with 


T := 


7 = T + Remainderft) 


((<m%\z-Z \ 2 + $\z-ri\%d x <p)' 


{(f> x , d\<j> x ) 


(49) 

(50) 

(51) 


(C) there exists a polynomial F(z,z) = 0(|^| 


such that the vector z 


satisfies the equation 

dt [\z\ 2 + F(z , z)\ = —T(z, z) + \z\Remainder[t) 
where r(z, z) is a positive quantity defined in (1191) . 


(52) 


The definition of R m ,n in d4lT) will be in Section [3 the proof of (l42l) will be 
in Section [9] C51) and (l50l) will be reformulated into Proposition llO.il (15211 will 
be proved in Section fill 


9 Proof of (02D 

Observe that in the equation for R in (1551) . the term on the right hand side, 
specifically JN(R,z ), contains terms quadratic in z and z. Hence for fixed 
z(t) £ C^, the equation for R(t) is forced by terms of order 0(\z(t)\ 2 ). 
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In what follows, we extract the quadratic in z,z part of R(t). Observe that 
the quadratic terms generated by the nonlinearity are of the form: 


^ JN m , n z m z n 

|m|+|n|=2 


2(f> x A 1 B 1 
-3 <j> x A\ - <j> x E 


(53) 


where A\ = afi, B\ = fi-r], and recall the definition of JN from (l28l) . Substitute 
this into the equation for R in (1551) and decompose R in the next results: 

Theorem 9.1. Define 

Rm,n ■- [L{ A) + iE( A) • (m - n) - 0 ]- x P c JN m>n , (54) 

and decompose R(t) as 

R = J2 R m,nZ m Z n + R (55) 

|m|+|n|=2 


Then the vector-function R(x,t) satisfies (1551) . 

The proof is the same to that in El, and skipping it will not affect under¬ 
standing the main part of this paper. Hence we choose to omit this part. 

To facilitate later discussions, we further decompose 

JN >2 = JN(R, z)- Y JN m , n z m z n . 

|m|+|n|=2 


We extract the third order terms of JN >2 '■ 

Y J N m,nZ m Z n = Y XRm,nZ m Z n +X 

|m| + |n|=3 |m| + |n|=2 

where X is a 2 x 2 matrix defined as 


A 2 

b 2 


and Ai, B\ and A 2 , 


( 2 <t> x B l 2 cf x A 1 

' ^ -6<?!> a Hi -20 a Hi 

B 2 are defined in (1551) . 


(Af + BDB, \ 
-{Af + B^Ax ) 

(56) 


(57) 


10 


Normal Form Transformon, Proofs of 
and 


In this section we present the proofs of equations (1551) and (15U1) , governing A and 
7 , crucial to controlling the large time behavior. The main result is Proposition 

nnn 
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This part is different from EHH], in that we have to define a new normal 
form transformation, some of whose parameters are defined as solutions to sys¬ 
tems of linear equations and their existence has to be justified. Moreover some 
small denominators will appear and we have to prove the numerators are also 
small. 

Now we present the idea. Central to our claim about the large time dynamics 
of NLS, is that the solution settles into an asymptotic solitary wave, (j} Xa ° , where 
A(t) —> Aoq. We achieve this by showing |A(f)| < eo(l + for some <5 > 0 

and small eo > 0. Since we expect the neutral mode amplitudes, z(t), to decay 
with a rate t~i , we require that there be no 0(\z(t)\ 2 ) on the right hand side 
of the equation (TTTTT 1 ) : 



The strategy is to choose the quadratic part of the polynomial a\(z,z) so as 
to eliminate all quadratic terms. There are two types of such terms, (1) the 
terms ZkZi and z k zi, and they are oscillatory with frequencies ~ —E k — Ei or 
~ E k + Ei, which stay away from zero. And the margins are large enough so 
that we can easily remove them by a normal form transformation, utilizing that 
ZkZi s=B dt _ i(E l +E ^ ZkZi. (2) The terms z k z m have frequencies ~ -E k + E m , 
which might be of small, or zero, frequencies. The key observation is that is 
that if the frequency of z k zi is small, then the coefficient is proportionally small! 
This allows us to define normal form transformation. The calculation is carried 
out below; see Lemma fit). 2 1 especially (l6T|) . 

Similarly we choose p n , q n , n = 1, 2, • • • , N, and 02 to remove most of the 
lower order terms in the equations for 7 , a n and /3 n . It turns out some terms can 
not be removed, for example the term \z n \ 2 in the equation for 7 , and \z n \ 2 z n 
terms in the equation for dtz. On the other hand, there terms either play a 
favorable role in our analysis, or will not affect it, namely it does not matter if 
7 is not convergent, since e* 7 ^ is only a phase factor. 

In defining p n , q n we have to solve system of linear equations, see e.g. (USD. 
The existence and uniqueness of solutions have to be addressed. 

In what follows we use the notations Nf™ n . N^ e n to stand for functions 
satisfying 



where, recall the definition of \m\ + \n\ = 2,3, from (l53l) and (l56l) . 

In what follows define the polynomials ai, 02 , p k and q k , k = 1,2, • • •, IV in 
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11221) by defining their coefficients: 


ai(z,z):= J2 ^] n Wz m z n , 

|m| + |n|=2,3 

02 {z,z)~ E A% n (\)z m z n , 

|m| + |n|=2,3, |m|^|n| 

Pk(z,z):= J2 PL%Wz m z n , (58) 

|m| + |n|=2,3 

q k (z,z):= Q%iWz m Z n , 

|ra| + |n|=2,3 


with m,n £ (Z + U {0})^. In what follows we use the notation 


m ■ E( A) = m kE k {\). 

k 


We start with defining Am] n - For \m\ = 2,3 


^m, 0 


« m o> A >- 


im ■ E( A) {(j) x 7 d\(t> x ) 


For |to| = 2, |n| = 1, 

4(!) — _1_ x 

m ’ n ' (0 A ,d A </> A )i(m-n)-F;(A) 

< / >A ) — ^ E T fc , n (r-r,^ A )], 

|fc| = |r| = |n| = 1, 
k + r = m 


(59) 


(60) 


here T mj „ is from the expansion of Y = T (z, z), defined in (ISTI) : 


T = E T rn,nZ m Z‘ 

|m| = |n| = l 


and the vector r is in (Z + U {0}) w . 
For \m\ = \n\ = 1, we define 


: = 


1 




{(j) x , d\(j) x ) i(m — n) ■ E(X)' m ’ n ’ 


(61) 

(62) 


Here in the first line it is possible that the denominator (m — n) • E(X) equals 
to zero or arbitrarily small, for example m = n, which might make the term 
ill-defined. In the second line we indicate that, by using Lemma 110.21 below, if 
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the denominator is small, then the numerator is proportionally small. Hence 
-^m|n &re always well defined, a similar calculation was in US- 

After defining various terms above, the other terms in Am] n are determined 
by the relations 

4^=4Sf°r|fc| + |/|=2,3. 

Now we define various terms in Amin, \m\ + \n\ = 2, 3, \m\ ^ \n\. 

For \m\ = 2,3, \n\ = 0, we define A^ 0 as 


~im ■ E(X)A%] 0 + := (0A> {N% 0 , d x cf> x ). 

For \m\ = 2, Ini = 1, 


(63) 


- i(m - n) ■ E(\)A%] n + A^ 

'' = {(t> x ,d x ^ N ™' nldA< ^> “ 2 


E 


T fc ,„(r-£,d A <^)], (64) 


| /c| = |n| = |r | = 1, k-\-r=m 


here the definition of T mjn and the convention on r are as in HMD- Especially 
in solving for A y m i n we need ( m _ n yE(\) t° uniformly bounded, indeed this is 
true by the facts Ek{ A) ~ eo — e k in (fT2l) and 2 ek < eg, see Condition (NL) in 
Section [3] 

( 2 ) 

The other terms in A, j are determined by the relations 


42“ 4ffcfor|fc| + |/|=2,3. 

Next we define coefficients Pmm and Qmin for the polynomials pu and qk, k = 
1 ,2, • • • ,N, see ([58]) . 

For |m| = 2, 3, |n| = 0, we define P^o and Qmo t° solutions to the linear 
equations 

-im ■ E(X)P. i fc) 0 - Ek(X)Qm] 0 := - Vk), 

-im ■ E(X)Q%] 0 + Ek(\)P^ 0 (65) 


Here the solutions exist and are unique because the corresponding 2x2 matrices 


-im-E(X) —Ek{\) 
Ek(\) -im-E(X) 


( 66 ) 


are uniformly invertible, or equivalently they determinants stay away from zero 
by a uniform margin. To see this, compute directly to obtain 

DetD = -[to • E(X)} 2 + El(X) = —[m • E(X) + E k (X)] [m ■ E(X) - E k (X)]. 

(67) 
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Next we relate the quantities on the right hand side to eo and ei, l = 1, 2, • • • , N, 
by (HU), namely 


for any l = 1,2, ■ ■ ■ , N, E t ( A) « e 0 - e*. 

Compute directly and use 2 ei < eo to see 

N 

m ■ E( A) - E k (X) ~ ^ m i ( e o - e i) - (eo - e fc ) 

i=i 

N 

= (|m| - l)e 0 + e fc - ^ 77 qe; ( 68 ) 

i=i 

is positive and stays way from zero for any m G U {0}} w and \m\ = 2,3. 
This together with that m ■ E(X) + E k ( A) is more positive implies the desired 
result: DetD stays away from zero. 

For |to| = 1 and \n\ = 2, we define Pm,n and to satisfy the equation 

„ p(fc) _/n(fc) 
m,n * 

+* [(»" ' V, Vk) ~ {r ' £, 6 )] 

|m|=|/e|=|r|=l,/e+r=n 

Ek — E(X) ■ (to — n) ’ 

(69) 

here the denominator E k — E{X) ■ (m — n ) stays away from zero by a uniform 
margin, by the same justification as in (1681) . and the definition of T m j j and the 
convention on r is as in (1601) . r • 77 := r k rj k and r ■ £ := Note that 

at this moment (1691) does not give unique solutions. This will be become clear 
in a moment. 

For \m\ = 2 and |n| = 1, we define 

iP^n ~ Q^n = 0. (70) 

After defining iPm}n — Qrh}n for (|m|, |n|) = (1, 2), (2, 1 ) above, it is not hard 

to see that these together with the relations Pm} n = Pn)m and Qm] n = Qn,m 
determine unique solutions for the linear equations. 

We continue to define Pm,h, Qm]n for |m| = \n\ = 1 , 

-i(m - n) • E{ A) - E k (X)Q™ n := - (N% n , Vk ), 

-*(to - n) ■ E( A) Q« n + E k {X)P^ n :=(N% n ,&). 

The solutions are well defined and unique since the matrix 

-i(m-n)-E{ A) -E k ( A) \ 

E k ( A) -i(m - n) • 2?(A) y 


(71) 

(72) 
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is uniformly invertible by the same arguments in showing the invertibility of the 
matrix in (16611 . 

We complete defining all the relevant terms by requiring that 


p(k) — p ( k ) (j(k) — p. 
r m.n r n,rm qr m n ■— q? 


(k) 

n,m* 


By now we have finished defining the polynomials ai, a- 2 , Pk- ( lk ; k = 
!, 2,- - - ,iV. 

Next we study the equation for i*.. By the definitions of coefficients Pm,L, 
we removed the following terms from the inequation: z m z n if |m| + |?r| =2,3 
and (|m|, |n|) (2,1). The result is: 

Proposition 10.1. Define the polynomials ai(z,z), a 2 (z,z), p n (z,z), q n (z,z ) 
as above. Then, (HE)-® holds and moreover for k = 1, 2, • • • ,7V 


+ iE k (X)z k 



(73) 


1, 


AT 


+ 2 T E^ 

m—1 


-iVm 


Vk 

i£,k 


+ Remainder (t). 


Proof. Recall the convention that Remainder represents any quantity satisfying 

< \z(tt + \\{x)- v Rmi + ii^wiil + m i n(*)-^(t)ii 2 . (74) 


We start with casting the A— and 7 —eqns in 


into a matrix form 


A 

7-T 


= fl + Remainder 


[Id + M(z, R,p, <?)] 

where, the vector fl is defined as 

( ImN ,(f X ) + | T ( {z-z)-r],(l) x ) }-d t a 1 


(75) 


n := 


pr§^[{ReN,d X (l) X ) - ( {z + z) ■ ^d X (f x )] - T - d t a 2 +a 1 


(76) 


the term Remainder is produced by the term y ^ ^ ^ ^ 

is the 2x2 identity matrix, M(z, R,p,q) is a matrix depending on z,R,p and 
q and satisfying the estimate 


|| M(z, R,p , q)|| = 0(\z\) + Remainder. 


(77) 


The smallness of the matrix M makes [Jd+M] uniformly bounded, hence 
by (1751) 


|A|, |7 — T| < \n\ + Remainder. 


(78) 
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Next we estimate f2, and start with casting it into a convenient form. 

The purpose of defining cii and <22 in (l58l) is to remove the lower order terms, 
in z and z, from ( ImN , (j) X ) — |T( (z — z) ■ rj, </> A ) and ( ReN , d X (j) X ) + |T( (z + 
z) ■ ^,d X (j) X ) to get 


n — D\ + D2 


(79) 


with 


D 1 := 


(<^ A , d\(j) x ) 


( —(ImN— IV, 

|ra| + |n|=2,3 

(ReN — J2 N % 


Im z m z n ,<t) x ) \ 


Re m —n 


Z Z 


\d x R) 


|m| + |n|=2,3 


and 


D 2—~ E 

|ra| + |n|=2,3 

It is easy to see that 


d t (A ( m}nZ m z n ) + iE( A) • (to - n)Am} n z m z r ‘ 
d t (Am] n z m z n ) + iE(X) ■ (to - n)A^] n z rn z ri 


D\ = Remainder. 


(80) 


To control D 2 we use a preliminary estimate from the z n — equation in (1351) 


d t z n + iE n (X)z n = 0(\z \ 2 ) + Remainder (81) 


to obtain 


D 2 


a E 

|ra| + |n|—2,3 


d x A&} n z m z n 
d x A%] n z m z n 


=0(\z\ 3 ) + Remainder 


+ 0(\z\ 3 ) + Remainder 


here the term 0(\z\ 3 ) is from the term 0(|z| 2 ) in dSTT) . 
Collect the estimates above to obtain 


(82) 


A, 7 — T =0(\z\ 3 ) + Remainder. (83) 

These estimates are still worse than the desired m-m- The reason is that 
their derivations depend on the non-optimal (1511) . Next we improve it using 


Choose p n and q n as in (1651) to remove the following lower order terms: z m z n 
satisfying \m\ + |n| = 2,3 and (|to|, |n|) ^ (2,1), to obtain 


dtZ n T iE n (X)z n 



+ Dz(n) + Remainder 


iz ■ rj 


Vn 


(84) 
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where £> 3 ( 71 ) is defined as 

A»(") “ - E M p k!l^) + *( fc - 0 ' 

|k|+|J|=2,3 

-i E + 

|k|+|i|=2,3 

= -A E Q&V* 1 

|fe| + |J|=2,3 

- E [ p i? - Qm } ] [ 5 *(— ) + *( fc - 0 • E{\)z k z l ] 

|k|+|J|=2,3 

=rr + r 2 (85) 

and Ti and r 2 naturally defined. 

For Ti the estimate for A in (1551) implies that 

Ti = Remainder. ( 86 ) 

For r 2 , the preliminary estimate in (1811) implies 

r 2 = 0(\z\ 3 ) + Remainder. 

This, in turn, implies an estimate better than (1811) 

dtz n + iE n (X)z n = 0(\z\ 3 ) + Remainder. (87) 

Sine the estimates derived for T 2 depends on non-optimal (1811) , this optimal one 
enables us to find 

r 2 = Remainder. (88) 

This together with Ti = Remainder in (1861) implies 

£>3 (n) = Ti + r 2 = Remainder. 

Put this into the dtz n -eqn in (1551) to obtain the desired estimate GS). 

(157)) also helps us to improve the estimate (1551) for D 2 

D 2 = Remainder. (89) 

This together with (1501) . (1751) implies for Q in (1751) 

fl = D i + D 2 = Remainder. 

Put this into (1751) to obtain the desired estimates (1551) . (15U1) for A and 7 — T. 

□ 

The following result has been applied in (l 6 lT) to show that the numerator is 
proportional to the denominator. Similar result can be found in HU- 
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Lemma 10.2. For \m\ = \n\ = 1, and m, n G (Z + U {0})^ 

(N^ n , (j> x ) = j.(n - m) ■ E( A) [(to • £, n ■ £} + (to • r/, n ■ rf)]. (90) 

Proof. We start with deriving an expression for \m\ = \n\ = 1. 

The explicit form of E JN rn ^ n z m z n in (l53l) implies that 

|m| + |n|=2 


E N L m n z m z n =2 <t> X A 1 B 1 

| ra | + | n |=2 

^ N N N N 

= ^ Zn£,n + ^ ^ Znfjn) ( ^ ^ Zm'Fjm ^ ^ ZmJin 

n—1 n= 1 m=l m =1 


Take the relevant terms to obtain 


N ^n = ~ ZruVn), 


(91) 


here we used the notation 


£n = n - Vn = n - rj, for n G (Z + U {0})^ and |n| = 1. (92) 

Hence the left hand side of (1901) takes a new form, 

(N%n, <f X ) I - UVn) 

=1[([L + (A) - L_(A)]C„,77m) - <[L+(A) - L_(X)\U,Vn)\, 

where we used the fact that L + (A) — L_(A) = 2(</> A ) 2 , see (flQ|) . Use the facts 
L + (A) and L_(A) are self-adjoint, and 

Z/_(A)7y n — .£/ n (A)£ n , and Z/_)_(A)^ n — £' n (A)?y n , n — 1,2, • • •, ./V, 

in urn and hence the desired result, recall the notations in (1^1) . 

(NL m n ,<j> X ) A) - £WA))[<£ m , Q + (Vm, Vn)] 

=-t(n — m) • -E(A) [(m • £, n ■ £) + (m ■ rj, n- if}] (93) 

1 / 

□ 
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11 Proof of the Normal Form Equation (1521) 


We expand the first two terms on the right hand side of the equations for z n in 
Gl to obtain 


- / E JN ™>‘ 

\ |ra|=2,|n|=l 

1=1 


m-n / Vk 

* ’ V -i£k 


1 


JV 


t E 

m= 1 


-iVn 


Vk 

’ V i£k 


(94) 


where, recall the definitions of JN m , n , \m\ + \n\ = 3, from (155jl . 0i(fc) is defined 
as 

e,(k) := - (xPfc.c ( ) ) = - ^ E W)* ( -i & ) ). 

(95) 

where, recall the definition of 2 x 2 matrix X in (1571) and we divide it into two 
terms X = X\ + X\ with Xi defined as 


, Z-T], (f) X z ■ £ 

1 ‘ ' —3 (j> x z ■ £, icj> x z ■ r] J ’ 


(96) 


0 2 (fc) := - ( Xj. 


V k 

k 


- x 1 


/ iV \ 

E E 

n=l |m| = |Z| = l 
JV 

E E 

V n=l|m[-|«l=l / 

/ N \ 

|m|=2 n=1 

\ |m|—2 n= 1 / 


Vk 

' V 


, (97) 


«•« : — s <«*■«’- «* * ( -^-’e )•( ^ 

+ i(( l ,xp + „.,n(r-).(* & 


*<*>=-5 T ((.T’M2i 
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where, recall the definition of real function T from (f50l) . 

The result is the following 

N 

Proposition 11.1. Re Zk®i(k) can be decomposed into three terms, 
k =1 


N 

Re ^2 z k @i(k) = - C(e 0 - A) T{z, z ) 

fc=l 

+ (e 0 — A) ^ C m , n (A) E(X)-(m-n) z m z n (98) 

|m| = |n|=2 

+ |e 0 — A| 2 0(M 4 ) 

where C is a positive constant, T(z,z) is the positive term defined in Fermi- 
Golden-rule (HU). C m ,n( A) are uniformly bounded constants and recall that eo — 
A > 0, in Lemma r m 

N 

For Re 2&0 2 (h), there exist some constants -D m , n (A) such that 
k= 1 

N 

Re Y, ^© 2 (*) =(eo - A) ^ D ffl ,n(A) £(A) • (m - n) z™z n - (99) 

fc=l \m\ = \n\=2 

For n = 3,4, 5, we have 


N 

Re^z k e n (k) = 0 . ( 100 ) 

fc=i 

The proposition will be proved in subsequent subsections. 

Next we prove the desired result (15^1) . 

Proof of (15^1) By the estimates above we have that 

d t \z\ 2 = - 2C(e 0 - A) T (z, z) + (e 0 - A) ]T C m , n E{ A) • (m - n)z m z n 

\m\ = \n\=2 

+ 0((e o — A) 2 | 2 | 4 ) + \z\Remainer, (101) 

with C mtn being some constant. Here we can take E Cm,n T/(A) • (jfl 

\m\ = \n\=2 

n)z rn z n to be real since it is in the equation for real parameter \z\ 2 and CT(z , z) 
is real. This implies 

^ C m , n E(\)-(rn-n)z m z n = ^ E(\) • (m - n)z m z n , 

|m| —|n|=2 |m| —|n|—2 
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and hence forces 


Cm,n — Cn,ra* 


( 102 ) 


By observing that 

E{ A) • (to — n)z m z n = idtz m z n + \z\Remainder , 
the fact that (eo—A) ^ iC™ „z m z n is real implied bv (11021) . we can define a 

|m| = |n|=2 


new nonnegative parameter z satisfying z 2 = \z\ 2 
such that 


(e 0 -A) Y, iC m,nZ m z n 

\m\ = \n\=2 


d t z 2 = —2C(eo — A) T{z,z) + 0((e o — A) 2 |z| 4 ) + \z\Remainer 
which is the desired estimate (1521) . 


(103) 

□ 


Next we prove Proposition 111.11 In the proof the following results, from 
m, will be used. Recall that the function <f> is the ground state for —A + V 
with eigenvalue — eo, the functions £(. m , k = 1, 2, • ■ ■ , N, are neutral modes with 
eigenvalues — e*,. 

Lemma 11.2. There exist constants Co, Ci £ R such that in the space (x)~ 4 TL 2 

4> x =Co(eo — A) 2 </> + C(|eo — A| 2 ), 
d x ^ x =Ci(eo - X)~^ + O(\e 0 - A|*), (104) 

---< 1. 

(</> a ,5a^ a ) ~ 

For the neutral modes we have 

m 4 (Vm - &\\ w 2 , ||<*> 4 (6n - &\\n^ IK*> 4 (£ m - = O(|eo - A|). 

(105) 

Recall Pis the orthogonal project onto the essential spectrum of — A + V, and 
P x of (fTT>)) is the Riesz projection onto the essential spectrum of L{ A) 

Pc = Pc n ( l J ) +O(|eo-A|). (106) 


11.1 Proof of (198|1 

We start from the definition of 0i(fc) in (f)5l) . Compute directly to obtain 

(IT)* ( Vk \ = ( ~ i( t )Xz ' ? ? dk + 3 i(j) x z ■ £ £ k \ 

1 V ) \ 4 X z ■ £, Pk + 4> Xz ■ d £fc ) ' 
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Hence 


— i(j> x (z ■ r]) 2 + 3 i(f> x (z ■ £) 2 
2(j) x z ■ £ z- r] 


i2e^f fc 0 1 (fc) = -/ ^ R m ,oz m , 

k \ |m|—2 

We extract its main part by define 


D := Re ^ < z m [(-A + V + A) J + iE( A) • m - 0] P c <j>(£ Hn y 

\m\—2 


(107) 


— l 

1 


<j){z ■ C Un ) 2 ' * 


>■ 

Here (£ hn ) m i s defined as, for m = (mi, m 2 , • • • , mjv) € (Z + U {0})^, 
^linym _ # . ^glinyrriN 


(108) 


The result is 

Lemma 11.3. 


N 

Re^2z k Qi(k) = -C(e 0 - A) D + 0(|e o - A| 2 |*| 4 ). (109) 

fc=i 

The lemma will be proved in Part 111.1.11 below. 

Next we study the term D. To facilitate our estimate we diagonalize the 
matrix operator in (11081) . Define a unitary matrix U by 


U := 


1 

72 



( 110 ) 


then we have that 

J = iUa 3 U*, 

with <73 being the third Pauli matrix. 

Inset the identity UU* = U*U = Id into appropriate places in the inner 
product of D to obtain, recall the convention that (/, g) = J f(x)g(x ) dx, 

D= 2 Re ( [*(- A + V + A - E{ A) • m) + 0] _1 P c <^ in ) m , 0(£ ,in )"> z m z n 

\m\-\n\=2 

=2 Im{[-A + V + X- E( A) • m - iO] _1 P c ^(^ iri ) m , 0(f I<n ) n ) 2 m ^ n . 

\m\ = \n\=2 

( 111 ) 

To cast the expression into a convenient form, we use the following two simple 
facts, for any functions /, g and real constant h, Im(f , g) = ^[(/, g) — {g, /)] 
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and (/, (—A + V — h — iO) x g) = ((—A + V — h + zO) 1 f, g). Compute directly 
to obtain 

D=- ^ Hn ) n ) z m z n . ( 112 ) 

M=M=2 

with L(m , n) being a linear operator defined as 

L(m, n) := [—A + 1/ + A — -E(A) • m — zO]^ 1 — [—A + V + A — E( A) • n + zO] -1 

In studying ( 11121 ) . the main tool is a well known fact that, see e.g. [IB] , for 
any constant h > 0, 

(y [[-A + V - h 2 - zO]- 1 - [-A + V - h 2 + zO]" 1 ] P c f, g) 

=Ch f f{ha) 'g(ha) da (113) 

J s 2 

where C is a positive constant, and the complex function / is defined as, 

f(k) ■.= [ f(x) e(x, k) dx, ( 114 ) 

Jr 3 

hence f(ha) is / restricted to the sphere |fc| = ha, a £ S 2 . Here the complex 
function e : K 3 x ]R 3 —> C is defined as 

e(x, k) := [1 + (—A - \k\ 2 - zO) _1 V(a;)] _1 e“' fe . 

We continue to study ( 11121 ) . For the easiest cases E(X) ■ m = E( A) • n, apply 
(llldD directly to obtain 

-(L(m,n)P c ^ Hn ) m , ^ Hn ) n ) = CC m [ <H^) m {C m a) <f{ip) n {C m a) da 

1 J s 2 

(115) 


with C m £ C defined as 

C m := V E (X) ‘ m ~ A, 

here E( A) ■ m — A is positive by the conditions that 2e*, < eg and -E&(A) ~ 
eg — &ki k = 1, 2, • • • , N, and A ~ eg. Recall that m £ (Z + U {0})^ and \m\ = 2. 
For the cases E( A) ■ m ^ E( A) • n, we claim, for some constant C m , n 

-(L(m,n)P c <t>(Z lin ) m , <f>(t lin ) n ) =cdck f <j>(t^ m (C m a) <H^ n (C n a) da 
1 Js 2 

+ C m>n E{\) ■ (m — n). ( 116 ) 

If the claim holds, this together with (11151) . Lemma Til. 31 and the fact Ek{ A) = 
eo — e-k + 0{e o — A) by (fT^l) . implies the desired result (1M1) . 
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What is left is to prove the claim (111611 . We start with decomposing the left 
hand side into two parts 


-(L(m,n)P c <i>(Z Hn ) m , <t>£ in ) n ) =A + B 

l 


( 117 ) 


with 


A :=-(L(m,m)P c <^r, </>(? m ) n ), 

l 

B :=-{[L(m,n) - L(m,m)\P c ^ Hn r, ^ Kn ) n >- 


By (gH , it is easy to see that 


A=CC m / <f>{t, lin ) m (C m a) <t>(£, lin ) n {C m a) dc t 

J s 2 

1 1 "ti 


( 118 ) 


=CC&C% / 4>{£, Un ) m {C rn a) (t>{e in ) n {C n a) da + C m , n E{ A) • (m - n), 

Js 2 

where C mtn is a constant, and in the second step we use that the functions 
(j>(^ hn )m an d the scalar C m depend smoothly on E(X) ■ m and E(X) ■ n. 

For B 1 it is easy to see that for some constant D m ^ n , 

B = D m>n E{ A) • (to - n). 

Collecting the estimates above, we prove the claim (11161) . 

Hence the proof is complete. 


11.1.1 Proof of Lemma 111.31 

We rewrite the expression in (11071) as 


N 


ReJ2zk®i(k) = - ( R ™,oz m , A 

\| m |=2 / 


k =1 

with the vector function A defined as 


A , = ( -i(t> (z ■ vf + 3 ■ £)" 

2 (j) x z ■ £ z ■ rj 


Apply the estimates of </> A , £ and r) in (11041) and (11061) to obtain 


A = C 1 (e 0 -\)b(z-t; lin ) 2 ' * 


+ O((e 0 -A)^|^| 2 ), 


( 119 ) 


( 120 ) 


here the expansion is in the space (x) 4 £ 2 for some eo > 0, C\ > 0 is a constant. 
Now we turn to R m $, which is defined as 

Rm, 0 = [£(A) + iE{ A) • TO - 0] _1 P A 
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and for JN m fi we have, from (15^1) 


E JN mfi z m 

\m\=2 


\ ' ( ImNmfi \ m 

v V -ReN m: o J ~ 

| m |—2 v ' 


1 

4 


—2 i(f> x z ■ £ z ■ r] 

-3 cj) x (z ■ £) 2 + 4> x (z ■ v) 2 


( 121 ) 


Now we extract the main part of P? m0 and JN m 0 by applying the estimates 
of cj) x , £ and ij and P c in (11041) and (11061) . and use that 

[L{ A) + iE{ A) • m - 0] _1 = [(-A + V + A) J + iE{ A) • m - 0] _1 + O(e 0 - A) 

to obtain 


R m ,o=C 2 (e 0 -A)5 z m [(-A + F + A)J + ? £(A)-m-0] 1 P^(f T1 ) m ^ ^ j 

+ O((eo-A)i)^ m (122) 

where C 2 > 0 is a constant, the expansion is in the space (x) 4 C 2 . 

Put this and (11201) into (11191) to obtain the desired result. 

□ 


11.2 Proof of (1991) 


To illustrate the ideas we consider part of it, namely 

N 


/ N \ 


y^z fc 9 2 (fc) 



V 


|m|—2 n= 1 

N 


E ^ ro E^,U 
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The other part is similar, hence omitted. Compute directly to obtain 


y^z fc 9 2 (fc) 


/ N \ 

‘ E * m E p £o& 


|m|=2 n=l 
N 


E * m E<U 




2: • 77 
—iz • £ 




\ |to|= 2 n=l 

/ N \ 

E * m E p S^ / J2 ReN ^~ \ 


=4 


1*1=2 

E ImNi fi z l 

V ife 


|m|=2 n=l 

JV - * 

E * to Eq£U. 

\ |m|—2 n=l / 

N 

= ^ 4j E E E [ P S*”* &*> i?e^E)+Q^ m (^n, ZmiViE)] (123) 

n=l | m |=2 |i|=2 


where, in the second last step we used the definition of Re Ni x, and ImN^ 0 , 
|Z| = 2, in (11211) . Next, relate the definitions of P©q an( i Qm ]0 i n (E3> to 
{rj n , ImNmfi) and (£ n , ReN mi 0 ), and then take the real part for Ylk^k® 2 (h) 
to obtain 


N _ _ 

Re E ^© 2 (fc) =4 E E Pe [I P S^S - Q %0 P i, o } ] 1 ■ E W z m z l ] 

k n= 1 |m| = |i|=2 

AT _ _ 

=2 e E Re p /,o } ] / • ^(A) * m i'] 

n=l |m| = |i|=2 

A/ - _ _ 

+ 2 E E Re iip^oQifi - Q%o p ifi] 1 ■ p w 

n =1 |m| = |Z|—2 

N _ _ 

=2 e E Re [\. p SoQ { $ - q2> p $] (i - m ) • p (A) 

n=l |m| = |i|=2 

= E (* - m ) ■ E W zm z l (124) 

M=l*l=2 

where, in the second step we used a simple fact that Re a = Re a, for any 
parameter a, for the second term in the second step, we interchange the indices 
m and l to obtain the third step, and in the last step C mt i are constants naturally 
defined. 

The proof is complete by observing that (I124D is the desired result. 
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11.3 Proof of (jinU 


It is easy to see that for j 
taking the real part, 


N 

= 3,4, ^^ZkQj(k) are purely imaginary. Hence by 
fc =1 


N 


Re ZkQj(k) = 0. 


(125) 


fc=l 


Next we turn to 0 5 . Use the definition of Xi and compute directly to obtain 

N / 


J2z k Q 5 (k) = -( Y R 

\ \m\ = \n\ = l 


Z rn z n 
m.n** & j 


k =1 


i(j> x \z ■ r)\ 2 + 3i(j) x \z ■ £| 2 
(j) X [z ■ £ z • 7] — z • t; • z • rj\ 


By noticing that R m ^ n z rn z n is real by definition, and that the other 

M=H=i 

vector function is purely imaginary, we have that the inner product is purely 
imaginary, hence 


N 

Re Y, = 0. 

fc=l 


(126) 


Collect the estimates above to complete the proof. 


12 Proof of the Main Theorem 16.1 


The proof is almost identical to the part in mi. hence we only sketch it. 

We begin by introducing a family of space-time norms, Z(T), 1Zj{T). for 
measuring the decay of the z(t) and R(t) for 0 < t < T, with T arbitrary 
and large. We then prove that this family of norms satisfy a set of coupled 
inequalities, from which we can infer the desired large time asymptotic behavior. 
Define 

r o :=|*(0)r\ (127) 

where, recall that |z(0)| is the initial amount of mass of neutral modes, see 
Theorem 16.11 

Now we define the controlling functions: 


Z(T) := max(T 0 + t) 2 |z(i)|, 
R. 2 (,T) := max(T 0 + £)||.R(t)|| 00 , 


U 4 {T) :=max||JS(t)|| W 3, 


Ki(T) := max(T 0 +t)\\(x)-'R(t)\\ n 3, 

1UT) :=max(T| +t)i||(x)-^(t)|| 2 , 

Ks(r) :=l gg ioIm+() l|jj(t)ll3 ~ 

(128) 
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To estimate 7 Zk, k = 1,2, 3,5, or to control the dispersion R, we use the 
propagator estimate, namely for any function g £ C 1 , 

l|e ti(A) P c A «?IU<t-^||5lli- (129) 

To estimate the decay of \z\, we dehne a positive constant q by 

q 2 = \z\ 2 + F(z,z) 

where F(z, z) is real and of order \z\ 4 , as stated in Statement C of Theorem l8.ll 
By Statement C of Theorem 18.11 

^M 2 < ~C\q\ A H-. (130) 

If the omitted part is small, then we have \q\ 2 = \z\ 2 < (1 +t) _1 . 

The details are identical to the corresponding parts in mi, and the proof is 
tedious. Hence we omit the detail. The results are: 

Proposition 12.1. 

Ki <T 0 ||(a;} ! 'i?(0)|| W 2 + 72 4 72 2 + Z 2 + T~ h [Z 3 + Z72i + 1Z 2 + 1Z 2 ]. 
n 2 <Tb||i?(0)||i + T o ||i?( 0)|| M 2 + z 2 + tz\tz 2 + T ~ 1 [Z 3 + zn x + n 2 }. 

1Z 3 <T 0 ||(a;) y i?(0)|| 2 +t||^| 2 (0) + T~* (Z 3 + Z1Z 3 + Z1Z i +K 3 + T 2 72 4 ). 
n\ < 11 ^( 0 ) 11^2 +T~ 1 [n 2 + z 2 Th + z 2 iz 2 + izliz 2 2 }. 

K 5 <T 0 ||i?(0)|| 1 + t 0 ||£( 0)||„2 + z 2 + T~--\n 2 n 2 + z 4 + zn 3 + n\ + n 2 2 \. 

Z{T ) <1 + 4r Z(T) ( Z(T) + Hl(T) + 7 Z 2 (T) + Z(T)1Z 3 (T) ) 

'T'b 

1 0 

12.1 Proof of the Main Theorem 16.11 

Our goal is to prove that the functions Z and 7 Zk, k = 1, 2, • • • ,5, are uniformly 
bounded, using the estimates in Proposition 1 12. II 

4 

Define M(T) := ^ 7 Z n (T) and 

n= 1 

5 := T o (|| 7?(0)||„2 + ||(x)^(0)|| 2 + \\m\\i) + ||^(0)||«2 +T<f N 2 , (131) 

where, recall the definition of To after (11281) . By the conditions on the initial 
condition HID we have that 72-4(0) is small, M{ 0), Z(0) and S are bounded. 
Then by the estimates in Proposition 112.11 we obtain 

Af(r)+ Z(T)</i(S), and 7*4 «1, (132) 
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where /r is a bounded function if S bounded. This together with the definitions 
of IZk, A: = 1,2,--- ,5, implies that 

II (x)~ v Rh, ||i?||oo <c(T 0 + f)-\ \z(t)\<c(T 0 +t)-i (133) 

which is part of Statements (A) and (C) in Theorem 16.11 The rest of (A), 
Equation (1231) . is proved in (1521) . 

□ 
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